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Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

inches 

2.54 

centimeters 

miles  (U.S.  statute) 

1.609 

kilometers 

pounds  (force) 

4.44822 

newtons 

pounds  (force)  per  square  inch 

6.894757 

kilopascais 

pounds  (mass)  per  cubic  foot 

16.011846 

kilograms  per  cubic  meter 

square  inches 

6.4516 

square  centimeters 

tons 

907.184 

kilograms 
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Introduction 


Background 

There  are  many  factors  used  to  determine  vehicle  trafiScability.  Some  factors, 
such  as  soil  strength,  may  even  require  field  measurement  or  be  inferred  based  on 
climatic  conditions.  The  NATO  Reference  Mobility  Model  (NRMM)  developed 
at  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Geotech¬ 
nical  and  Structures  Laboratoiy  (GSL),  Vicksburg  site,  uses  the  data  fi-om  over 
100  factors  that  describe  a  vehicle  off-road  terrain  unit,  road  unit,  or  linear  feature 
to  predict  potential  maximum  safe  vehicular  speeds.  These  factors  include  smface 
composition,  geometry,  vegetation,  and  linear  geometry  (Tumage  and  Smith 
1983,  Harrell  1989).  In  order  to  characterize  tactical  logistics  situations,  certain 
U.S.  Army  Training  and  Doctrine  Command  (TRADOC)  operational  scenario 
categories  are  used  with  NRMM,  such  as  tactical  high  (e.g.  >  30  kph)  or  tactical 
low  (<  30  kph)  to  predict  vehicular  speeds.  Such  analyses  have  been  used  previ¬ 
ously  by  the  ERDC,  TRADOC,  and  the  acquisition  community  (U.S.  Army  Tank- 
automotive  and  Armaments  Command  (TACOM))  to  evaluate  the  performance  of 
vehicles  in  given  scenarios.  Some  of  the  goals  of  the  AT40  6.2  Research 
Program,  “Lines  of  Communication,  Repair  of  Roadways”  are  to:  (a)  develop 
methodologies  and  analytic  software  for  determining  appropriate  convoy  speeds, 
tire  pressure,  and  load  capacities  on  roadway  components,  and  (b)  develop 
methodology  to  calculate  throughput  as  a  fimction  of  visibiility  and  weather 
conditions.  One  of  the  critical  parameters  that  governs  speed  predictions  for  a 
given  vehicle  is  its  ride  response  over  areas  of  terrain  and  road  surface  roughness. 
During  force  projection  operations  into  underdeveloped  or  undeveloped  countries 
with  poor  transportation  networks,  logistics  throughput  and  maneuver  operations 
will  be  greatly  affected  by  speed  performance  and  driver  ride  conditions,  especi¬ 
ally  as  roads  deteriorate  with  traffic  and  changing  weather  conditions.  Thus, 
given  a  tactical  situation,  it  is  important  to  develop  a  means  of  understanding  how 
to  rapidly  assess  the  infiuence  of  deteriorating  road  surface  conditions  on  ride 
response  of  vehicles  that  are  used  by  the  force  and  evaluate  this  effect  in  the 
tactical  planning  process. 
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Purpose 

The  purpose  of  this  study  was  to  develop  and  recommend  a  methodology  for 
evaluating  tactical  support  road  situations  which  will  assist  in  determining  appro¬ 
priate  convoy  speeds,  evaluate  the  use  of  bypass  routes  around  lost  critical  nodes 
in  the  infrastructure,  and  allow  the  computation  of  road  load  capacities  for  deter¬ 
mining  the  logistics  vehicle  requirements  for  various  tactical  scenarios.  Thus,  the 
impacts  on  throughput  and  maneuver  can  be  evaluated  in  quantifiable  terms. 


Scope 

Sixteen  representative  field  experiments  with  a  prototype  Logistics  Support 
Vehicle  (LSV)  were  conducted  at  two  local  test  sites.  Mound,  LA,  and 
LeToumeau,  MS.  At  the  Mound  site,  the  test  vehicle’s  accelerometers  were 
calibrated,  and  two  sections  of  road,  one  asphalt  and  one  gravel,  were  used  to 
evaluate  the  technology  under  development.  At  the  LeToumeau  site,  four  rock- 
covered  ride  courses  were  used  to  collect  accelerometer  and  ride  meter  data  over 
courses  with  known  surface  roughness.  This  report  is  part  of  a  series  of  technical 
reports  plaimed  to  develop  and  describe  the  methodology  for  using  the  LSV.  The 
first  report  described  the  full  characterization  of  the  vehicle  dynamic  response  of 
the  LSV,  and  this  report  will  describe  the  methodology  developed  for  using  the 
LSV  to  assess  the  impact  of  road  surface  conditions  in  the  supply  network  on 
logistics  throughput.  In  this  report,  an  irmovative  methodology  using  wavelet 
transforms  is  developed.  These  are  used  to  characterize  the  vehicle’s  vibrational 
response  to  ride  conditions  while  traversing  nondeformable  areas  of  known 
surface  roughness.  Because  the  mathematical  background  necessary  to 
understand  the  application  of  wavelets  to  vehicle  vibrational  response  is  relatively 
new,  a  short  review  of  wavelet  mathematics  is  also  included  in  this  report. 
Historically,  the  ERDC  has  detrended  rod  and  level  elevation  measurements. 
Instead  of  using  ride  performance  of  a  known  vehicle  to  determine  surface 
roughness,  the  root  mean  squared  statistic  has  been  used  to  characterize  surface 
roughness  from  the  detrended  elevation  values.  The  development  of  wavelet 
technology  as  an  alternative  means  to  characterize  the  surface  roughness  of  the 
ride  courses  through  use  of  the  LSV  in  a  similar  manner  is  also  described. 
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2  Test  Vehicle,  Experiments 
Conducted,  and  Data 
Collected 


Logistics  Support  Vehicle  (LSV) 

The  Logistics  Support  Vehicle  is  a  militarized  research  version  of  the  y4-ton 
Chevrolet  Suburban  4x4.  Similar  vehicles  were  purchased  by  the  military  under 
the  Commercial  Utility  Cargo  Vehicle  program  for  general  and  administrative  use 
by  ERDC-Vicksburg  personnel  to  measure  accelerations  at  each  wheel  and  the 
driver’s  position,  in  a  research  effort  to  determine  the  effects  of  deteriorating  road 
surface  conditions  on  vehicle  response.  Alternatively,  the  results  could  be  used  to 
estimate  road  surface  deterioration  and  road  life  as  a  function  of  the  traffic  speeds, 
wheeled  tire  pressure,  etc.  These  data  are  required  for  the  planning  of  convoyed 
military  vehicles  operating  on  roads  or  trails.  Long-term  research  could  lead  to 
procedures  for  characterizing  road  surface  changes  on  the  known  response  of  a 
specific  vehicle  used  as  a  calibrated  surrogate. 

The  LSV  is  a  four-wheel  drive,  diesel-powered,  4x4  carry-all  weighing 
1,635  lb  at  the  right  fi'ont,  1,570  lb  at  the  right  rear,  1,740  lb  at  the  left  front,  and 
1,685  lb  at  the  left  rear,  for  a  total  gross  vehicle  weight  of  6,630  lb.  The  vehicle 
instrumentation  consists  of  seven  accelerometers  located  on  the  vehicle  and 
labeled  as  follows:  Chi  =  Front  Left  Axle,  Ch2  =  Front  Left  Frame,  Ch3  =  Rear 
Left  Axle,  Ch4  =  Rear  Left  Frame,  Ch5  =  Front  Right  Axle,  Ch6  =  Front  Right 
Frame,  Ch7  =  Under  the  Driver’s  Seat. 


Field  Experiment 

During  March  through  April  2000,  tests  were  conducted  at  two  sites  in  the 
Vicksburg  area:  Mound,  LA,  and  LeToumeau,  MS.  The  road  sections  at  Mound, 
one  asphalt  and  one  gravel,  and  four  standard  test  courses  at  LeToumeau  were 
used  for  these  tests. 
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a.  Mound,  LA. 


( 1 )  Asphalt  road  to  airport. 

(2)  Bridge  detour  road,  gravel. 
b.  LeToumeau,  MS. 

(1)  Ride  Test  Course  1,  rms*  =  0.77. 

(2)  Ride  Test  Course  3,  rms  =  0.50  (gravel  road). 

(3)  Ride  Test  Course  4,  rms  =1.16. 

(4)  Ride  Test  Course  6,  rms  =1.81. 


Mound,  LA 

At  Mound,  the  vehicle  was  operated  at  a  constant  speed  of  30  mph  over  the 
asphalt  road  and  12  mph  over  the  gravel  road,  while  vibrational  data  were  col¬ 
lected  at  each  accelerated  location.  These  data,  along  with  baseline  drop  tests 
conducted  with  the  LSV  at  ERDC,  Vicksburg,  to  be  reported  in  another  report  of 
this  series,  will  be  used  for  three  research  purposes:  (1)  characterize  the  dynamic 
response  of  the  LSV  such  that  it  could  be  used  as  a  tool  for  rapidly  identifying 
surface  roughness  (rms)^  after  being  driven  over  a  course  or  road  section, 

(2)  develop  criteria  for  use  in  predicting  the  ride  performance  of  a  new  or  concept 
vehicle  as  a  function  of  known  vehicle  parameters  related  to  LSV,  and  (3)  use  the 
LSV  baseline  response  with  periodically  collected  data  of  a  known  section  of  road 
to  estimate  the  failure  rates  of  road  surfaces  as  a  fiinction  of  traffic  type  and  traffic 
intensity.  These  research  objectives  are  being  accomplished  under  the  AT40  6.2 
Research  Program,  ‘"Lines  of  Communication,  Repair  of  Roadways,”  which  is  a 
Science  and  Technology  Objective  (STO  IV.EN.  1998.03) 


LeToumeau,  MS 

At  LeToumeau,  the  vehicle  was  operated  at  a  variety  of  constant  speeds  while 
traversing  ride  courses  1,3,4,  and  6.  Sufficient  tests  were  conducted  to  define 
the  6-watt  absorbed  power  ( )  ride  course.  These  data  were  necessary  to 

characterize  the  baseline  ride  dynamics  of  the  LSV  over  curve  of  known  surface 
roughness  (rms).  Results  of  these  tests  are  shown  in  Table  1. 


^  Root  mean  squared  (rms)  =  y  ^ ,  (x^  -  where  Xi  is  the  measured  profile 

elevation  in  inches. 
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Data  Collected 


For  these  dynamics  response  tests  with  the  LSV,  eight  accelerometers  were 
used.  Each  was  mounted^onded  to  a  metal  mounting  block  of  aluminum  and 
secured  to  the  vehicle  at  the  desired  location.  These  accelerometers  are  direct 
current  (DC)  servo-type  (piezo-resistive),  which  provide  excellent  stability,  high 
signal-to-noise  ratio,  and  a  flat  frequency  response  in  the  band  region.  The  typical 
output  sensitivity  is  5  to  10  mV/g,  depending  on  the  selected  g  range.  They 
provide  a  frequency  response  two  orders  of  magnitude  higher  than  typical 
vehicular  acceleration  frequency  levels  and  are  normally  5-  or  10-g  range  devices. 
The  10-g  accelerometers  are  mounted  on  the  frame  with  one  mounted  on  the  floor 
beneath  or  directly  behind  the  driver’s  seat  of  the  LSV. 
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3  Background  to  Wavelet 
Analysis 


What  are  Transforms? 

Transforms  are  the  result  of  transforming  a  function  of  one  independent 
variable  to  that  of  another.  A  common  example  of  a  transform  of  f  (t)  is  the  cosine 
transform 

/(0  =  2]4cos(2;r^)  (1) 

ifc=0  ^ 

or 

/(O  =  cos(2;r©i0  (2) 

i=0 


where  ro,  =  klT  is  the  frequency.  The  cosine  transform  of  fif)  is  literally  the 
coefficients  in  the  series  Ak  =  A(coJ  which  are  given  by 


4 


_2 

T 


T 

J  f  {t)  cos{2Ka)^t)dt 
0 


(3) 


Thus,  the  function  in  the  time  domain  is  transformed  to  a  function  m  frequency 
domain.  The  cosine  transform  is  an  example  of  a  more  general  class  of  linear 
transforms  referred  to  as  Fourier  transforms.  Note  that  the  transform  is  a  linear 
combination  of  the  sequence  of  cosine  functions,  (a  property  of  the  transform), 
and  the  constants  Ak  (a  property  of  the  data).  The  cosine  functions  are  referred  to 
as  the  basis  functions  or  kernel  of  the  transform.  Other  choices  of  basis  functions 
give  different  transforms  having  different  properties. 

Consider  a  function  as  a  list,  or  time  record,  of  N  numbers  that  is  a  discrete 
sampling  of  a  continuous  function  in  time,  flf)  .  The  readings  are  equally  spaced 
over  a  total  time  of  T  =N  At  such  that  the  position,  ,  in  the  list  corresponds  to  a 
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time  equal  to  klsx .  In  the  discrete  transform,  the  list  of  N  samples  in  time  is 
transformed  to  a  list  of  N  numbers,  each  corresponding  to  a  discrete  frequency. 
Each  record  can  be  computed  from  the  summation 

fn=J^A  cos(2;rm^«A0  (4) 

k=0 

For  the  general  discrete  Fourier  transforms,  the  basis  functions  consist  of  both 
sines  and  cosines, 

N 

f„=^ [A^  cos(2jt(o^nAt)  +  (5) 

*=0 


which  is  generally  written  in  complex  form  as 

k=NI2 

fn= 

k=-NI2 

The  algorithm  for  computing  this  transform  efficiently  is  referred  to  as  the  Fast 
Fourier  Transform  (FFT). 


Why  Use  Transforms? 

The  transformed  data  carry  the  same  information  as  the  original  record, 
whereby  given  the  transform  the  original  data  can  be  retrieved.  However,  the 
transformed  data  are  sometimes  a  more  useful  representation.  For  example,  the 
resonant  frequency  of  the  vehicle  and  the  roughness  of  the  road  might  dominate 
the  record  of  a  road  test.  The  time-domain  record  might  not  show  the  superposed 
effects,  whereas  the  two  effects  would  be  obvious  in  the  frequency  domain.  Fil¬ 
tering  the  effect  of  roughness  from  the  record  amounts  to  removing  those  data 
from  the  frequency  domain  and  transforming  the  result  back  to  the  time  domain. 
Data  compression  techniques  make  use  of  the  dominance  of  particular  frequen¬ 
cies  in  a  record.  Compression  is  achieved  by  storing  or  transmitting  only  the  Fk 
that  are  significant;  decompression  then  is  the  transforming  back  to  the  time 
domain  to  obtain  an  approximate  record  that  is  now  missing  the  (presumably) 
irrelevant  data.  Unfortunately,  not  all  transforms  are  equally  efficient.  For 
example,  a  record  that  contains  only  one  nonzero  value  (a  spike)  is  most  effi¬ 
ciently  represented  in  the  original  time  domain.  In  the  frequency  domain  repre¬ 
sentation  of  a  spike,  each  element  of  the  Fourier  transform  has  the  same  value. 


What  are  Wavelets? 

Wavelets  are  mathematical  families  of  functional  atoms  or  functional  com¬ 
ponents.  They  can  be  generated  from  a  single  function,  *F,  by  dilations  and 
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translations  in  the  form:  =  \a\  -b)/  a);  where  the  parameters  a  and  b 

are  often  restricted  to  a  discrete  set.  In  this  form  they  provide  a  way  to  decompose 
an  arbitrary  function  similar  to  the  way  sines  and  cosines  are  used  in  the  mathe¬ 
matical  theory  of  Fourier  analysis.  Wavelet  Transforms  are  defined  and  used  in 
wavelet  methods  to  analyze  signals  in  a  way  analogous  to  Fourier  Transforms: 

For  a  and  b,  arbitrary  positive  or  negative  real  numbers,  the  discrete  wavelet 
transform  (DWT)  is 

(Tf)„  =  {¥„,„,/)  -  \dt  -  «^o)  *  /(O  (7) 


where  a  =  b  —  nb^a^ 

These  methods  provide  an  improved  way  to  analyze  and  represent  fiinctions  in  the 
spatial  form,  the  representations  of  which  are  changing  with  time.  The  book  by 
Krantz  (1999)  gives  a  good  overview  of  some  of  the  mathematical  topics  in  the 
theoiy  of  real  analysis  and  functional  analysis,  which  played  a  part  in  their 
invention  and  development. 

The  uses  of  wavelets  are  varied.  Fournier  (1995)  contains  an  introduction  to 
their  use  in  computer  graphics  and  image  compression,  and  Charles  Chu  (1997) 
explains  their  uses  to  improve  the  functional  representation  by  interpolation  of 
splines  of  data  for  numerical  partial  differential  equations.  The  papers  by 
Newland  and  Butler  (1999a,b)  use  wavelets  for  the  study  of  centrifuge 
experiments  related  to  earthquake  engineering.  The  reference  by  Torrence  and 
Compo  (1998)  provides  data  on  using  them  to  prove  the  statistical  significance  of 
El  Nino  related  sea  surface  temperatures. 

And,  corresponding  to  their  different  uses,  there  are  many  different  types  of 
wavelets:  continuous  wavelets  (Meyer  and  Ryan  1993),  discrete  wavelets  (Chu 
1997),  orthogonal  wavelets,  biorthogonal  and  compact  wavelets  (Daubechies 
1992).  The  Daubechies  scaling  function,  wavelet  function,  and  filter  coefficients 
are  shown  in  Figure  1 .  The  definitions  of  the  scaling  function  and  wavelet 
function  are  based  on  the  mathematical  algorithm  used  to  compute  the  wavelet 
decomposition  by  Chu  1997.  The  wavelet  filter  coefficients  are  the  multipliers 
used  in  the  wavelet  transform  to  extract  the  high-  and  low-frequency  components 
of  the  signal.  Chu  1997  has  a  good  reference  to  help  understand  how  wavelet 
function,  scaling  function  and  filter  coefficients  work  together  to  decompose, 
reconstruct,  and  approximate  the  signal. 

The  wavelet  decomposition  process  is  initiated  by  first  separating  a  signal 
into  two  parts,  an  approximation  part  A1  and  a  detailed  part  D1 .  This  is 
accomplished  by  using  the  coefficients  of  a  high-  and  low-pass  filter,  which  multi¬ 
ply  the  wavelet  and  scaling  function.  In  the  next  level  (step),  the  approximation 
A1  is  itself  separated  into  two  parts,  A2  and  D2.  This  process  decomposition  is 
repeated  until  the  final  level  of  wavelet  approximation  is  reached. 
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Scaling  function  phi 
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Wavelet  function  psi 
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Figure  1.  Example  of  wavelets 

This  report  will  primarily  consider  Daubechies  (1992)  wavelets  and  Meyer 
and  Ryan  (1993)  wavelets  as  defined  in  the  MATLAB  wavelet  toolbox  (Misiti 
et  al.  1996),  and  harmonic  wavelets  as  defined  by  Newland  (1993).  The  harmonic 
wavelets  were  picked  for  their  usefulness  in  analyzing  time  series  data  from 
vibrational  studies  such  as  those  reported  in  the  following:  “For  the  purposes  of 
signal  analysis,  a  wavelet  is  considered  to  be  a  bandpass  window  function  that 
stops  at  least  the  zero  frequency”  (Chu  1 997). 

Wavelet  packets  are  a  generalization  of  wavelets,  which  use  two  high  and  two 
low  bandpass  filters  instead  of  one  each  (Misiti  et  al.  1996,  Mallat  1998,  Meyer 
and  Ryan  1993,  Chu  1997).  Wavelet  packets  are  more  often  used  for  one¬ 
dimensional  (1-D)  signal  analysis  (as  opposed  to  two-dimensional  (2-D)  image 
compression),  because  they  do  a  better  job  of  representing  the  time-dependence  of 
the  signal.  They  use  a  family  of  two  sets  of  orthogonal  bases  to  decompose  the 
signal,  instead  of  one.  An  optimal  base  inside  of  the  family  can  be  chosen  that 
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minimizes  the  sum  of  squares  of  residuals  of  the  signal  values  and  the  wavelet 
approximation.  This  optimal  wavelet  base  is  often  represented  in  terms  of  a  par¬ 
ticular  wavelet  decomposition  tree.  The  wavelet  packet  decomposition  process 
begins  similarly  to  the  wavelet  decomposition  process.  First,  the  signal  is 
separated  into  two  parts,  an  approximation  part  A1  and  a  detailed  part  Dl.  This  is 
accomplished  by  the  coefficients  of  high-  and  low-pass  filters,  which  multiply  the 
wavelet  and  scaling  function.  Now,  in  contrast  to  the  wavelet  decomposition 
process,  in  the  next  level  (step)  of  the  packet  decomposition,  the  approximation 
A1  is  separated  into  two  parts,  A21  and  A22,  and  the  detailed  Dl  is  separated  into 
two  parts,  D21  and  D22.  This  is  a  composition,  as  mentioned  above,  by  using  two 
orthogonal  bases  for  the  decomposition.  One  base  corresponds  to  the  scaled  and 
translated  time/frequency  atomic  terms  generated  by  the  “mother”  wavelet 
function.  The  other  base  corresponds  to  the  scaled  and  translated  time/frequency 
atomic  terms  of  a  dual  “father”  wavelet  function.  The  process  is  repeated  until  the 
final  level  of  wavelet  approximation  is  reached. 

The  2-D  plot  of  the  square  of  the  signal  approximation  in  space  and  in  time  is 
called  a  Heisenberg  box  (Mallat  1998).  These  plots  can  be  examined  to  determine 
the  best  way  to  subdivide  the  spatial  and  temporal  axes  to  compute  the  wavelet 
bases. 

Figure  2  shows  a  wavelet  decomposition  tree  and  the  Heisenberg  box  for  the 
decomposition  of  a  power  spectral  density  function.  The  x-axis  of  the  signal 
(figure  in  the  upper  right  box)  is  frequency  and  the  y-axis  is  the  square  of  the 
amplitude  of  the  Fourier  transform  of  the  starting  signal.  The  x-axis  of  the  box 
(figure  in  box  at  lower  right)  is  time  and  the  y-axis  is  frequency.  The  number  of 
boxes  on  the  vertical  y-axis  corresponds  to  the  number  of  frequency  decomposi¬ 
tion  nodes  on  the  lowest  level  in  the  wavelet  decomposition  tree. 


Reasons  for  the  Use  of  Wavelets  in  Vibrational 
Analysis  Studies 

The  tutorial  by  Polikar  (1999)  gives  a  beginner’s  introduction  for  engineers  on 
the  use  of  wavelets  for  signal  analysis.  And,  the  more  advanced  reference  by 
Daubechies  (1992)  contains  most  of  the  pertinent  basic  theorems  in  mathematical 
analysis,  which  are  necessary  to  justify  the  use  of  wavelets  to  analyze  time  series 
data.  Fourier  transformations  separate  different  parts  of  the  data,  which  cor¬ 
respond  to  its  different  cyclic  components.  They  allow  an  experimenter  to  exam¬ 
ine  the  behavior  of  a  system  in  both  the  spatial  and  frequency  domains.  In  signal 
analysis,  a  series  of  related  data  observations  taken  over  a  period  of  time  is  called 
a  time  series.  If  the  statistical  parameters  that  characterize  different  parts  of  the 
series  of  data  vary  with  respect  to  time,  the  time  series  is  said  to  be  nonstationary. 
(Otherwise,  it  is  said  to  be  stationary.) 

There  are  several  reasons  why  it  is  advantageous  to  use  Wavelet  transforms 
instead  of  the  more  familiar  Fourier  Transform  to  analyze  time  series: 
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Figure  2.  Example  of  wavelet  decomposition 

“A  disadvantage  of  Fourier  analysis  is  that  the  frequency  information 
can  only  be  extracted  for  the  complete  duration  of  a  signal.  Since  the 
integral  in  the  Fourier  transform  extends  over  all  time,  the  information  it 
provides  arises  from  an  average  over  the  whole  length  of  the  signal.  If, 
at  some  point  in  the  lifetime  of  the  signal,  there  is  a  local  oscillation 
representing  a  particular  feature,  this  will  contribute  to  the  calculated 
Fourier  transform,  but  its  location  on  the  time  axis  will  be  lost. .  .By 
breaking  the  signal  into  shorter  and  shorter  time  segments,  using 
windowing  functions  and  dilation  equations,  wavelets  analyze  and 
extract  this  local  information. .  .Fourier  analysis  fails  to  capture  the  non¬ 
stationary  aspects  of  the  time  series,  but,  the  wavelet  approach  allows 
these  factors  to  be  characterized.”  (Newland,  1993  (p.  294)). 
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Harmonic  Wavelets 


Harmonic  wavelets  are  certain  combinations  of  Gabor  wavelets  introduced  in 
the  book  by  Newland  (1993).  Harmonic  wavelets  are  Gabor  wavelets;  thus,  they 
are  wavelet  packets  as  we  defined  them  previously.  Dennis  Gabor  introduced 
Gabor  wavelets  in  1946  in  a  study  of  speech  analysis  and  reconstruction.  They 
facilitate  the  representation  of  time  signals  as  a  superposition  of  elementary 
wavelets.  If  the  signal  is  expected  to  have  a  certain  temporal  structure  as  well  as 
certain  emitted  frequencies,  it  is  necessary  to  analyze  the  signal’s  properties  both 
with  respect  to  its  frequencies  and  the  time  when  they  were  emitted.  As  explained 
by  Walker  (1999)  and  Meyer  and  Ryan  (1993),  a  Gabor  wavelet  of  width 
parameter  w  and  frequency  parameter  V  has  a  wavelet  function. 

y/(x)  =  *e{-7t{xl M>Y)*e{i*2K  *v*{xIm>))  (8) 

The  independent  variable  in  this  equation  is  a  real  variable.  Harmonic  wave¬ 
let  packets  have  a  wavelet  function  of  a  complex  variable 

=  (e(-4;r  x)  -  e(i  2;r  x))  /(i  2;rx)  (9) 

So,  they  are  a  combination  of  two  complex  Gabor  wavelets  of  frequency  param¬ 
eters  4  and  2  and  each  with  width  parameter  1.  A  disadvantage  of  using  this 
wavelet  is  that  it  decays  on  the  x-axis  only  proportional  to  1/x.  Thus,  the  theo¬ 
rems  in  Meyer  (1993)  and  Daubechies  (1992)  concerning  r-regular  wavelets, 
where  r  >  1,  do  not  apply.  This  wavelet  function  does,  however,  generate  an 
orthogonal  basis,  which  can  be  used  for  multiresolution  analysis.'  And,  in  many 
cases,  its  simple  analytical  formula  is  an  advantage.  Also,  since  the  wavelet 
function  is  complex,  only  one-half  as  many  coefficients  need  to  be  computed  to 
determine  its  value  at  each  translation  and  dilation  index. 

Figure  3  shows  a  wavelet  function  and  the  windowing  transforms  for  the 
wavelet  and  scaling  functions  for  a  set  of  harmonic  wavelets,  which  will  be  used 
in  this  report. 

If  regularity  and  vanishing  moments  of  the  approximation  flmctions  are 
important,  Meyer  wavelets  can  be  used.  These  wavelets  are  analytic,  symmetric, 
with  a  compactly  supported^  Fourier  transform,  and  their  size  decays  faster  than 
any  polynomial.  They  are  shown  in  Figure  4.  The  paper  by  Abry,  Goncalves,  and 
Flondrin  (1995)  discusses  the  advantages  of  using  wavelets  with  r  vanishing 
moments  to  perform  power  spectral  estimation.  The  book  by  Kay  (1988)  is  a 
good  introduction  to  the  general  area  of  power  spectral  estimation  and  the 
problems  involved  in  constructing  good  statistical  sampling  algorithms  for  it. 


‘  These  theorems  guarantee  that  the  square  norm  of  Fourier  transform  of  wavelet  function  is 
infinitely  differentiable,  and  that  the  first  r  moments  of  the  wavelets  coefficients  are  zero.  This 
makes  them  desirable  to  use  as  approximations  and  interpolations. 

^  A  function  is  compactly  supported  if  the  part  of  its  domain  of  definition  where  it  is  nonzero 
contains  all  its  limit  points  and  is  a  bounded  set  with  respect  to  distance. 
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Figure  3.  Harmonic  waveiets 


Scaling  function  phi  Wavelet  function  psi 


Figure  4.  Meyer  wavelets 


Meyer  Wavelets 

Yamada  and  Sasaki  1998  describe  the  use  of  Meyer  wavelets  (Figure  4)  for 
the  analysis  of  atmospheric  wind,  turbulent  fluid,  and  seismic  acceleration  data. 
As  mentioned  previously,  the  wavelet  function  decreases  to  zero  when  x  increases 
to  infinity  faster  than  any  polynomial.  This  property  also  holds  for  the  derivatives 
of  this  wavelet  function.  These  wavelets  are  analytic,  symmetric,  and  orthogonal. 
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4  Wavelet  Analysis  for  the 
Ride  Experiments 


General  Background 

The  subject  of  random  vibrations  as  they  relate  to  the  dynamic  projection  of  a 
vibrating  vehicle  (mass,  stiffness,  damping,  etc.)  and  the  understanding  of  the 
statistical  characteristics  of  the  motion  generated  by  operating  over  on-/off-road 
terrains  will  be  very  important  as  the  Army  transforms  into  the  21st  century. 
Typically,  analysis  of  vehicle  motion  transforms  produces  vibration  data  averaged 
over  sections  of  roads  or  trails.  Local  position  of  particular  features  or  sections  of 
interest  is  lost  with  Fourier  analysis.  Wavelets  are  a  recent  development  in  signal 
analysis,  which  allows  decomposition  of  the  signal  into  a  series  of  other  generally 
orthogonal  functions  of  different  lengths  and  positions  with  respect  to  time. 
Features  producing  extreme  or  unusual  vibration  can  be  located  using  the  posi¬ 
tions  of  the  wavelets  into  which  the  signal  can  be  decomposed.  In  this  manner, 
vehicle  vibrational  data  collected  in  the  time  domain  may  mask  the  most  distinc¬ 
tive  information,  which  is  hidden  in  the  frequency  content  of  the  signal  test.  The 
FFT  algorithm  produces  frequency  domain  coefficients,  which  describe  the 
amplitude  and  phase  of  the  signal’s  harmonic  components.  The  disadvantages  of 
analyzing  only  time-domain  data  or  frequency  domain  FFT  data  are  defeated 
through  use  of  wavelet  analysis,  which  maintains  both  time  and  frequency 
resolution.  The  application  of  this  relatively  recent  technology  to  the  LSV  is 
discussed  in  this  report. 

In  this  study,  raw  data  collected  with  the  LSV  on  roads  at  Mound,  LA,  at 
100  samples/sec  per  channel  were  stored  in  units  of  g.'  The  interactive  software 
package  MATLAB,  with  a  wavelet  toolbox  written  by  Newland  (1993),  was  used 
to  evaluate  the  data  collected  at  Mound,  LA.  Using  discrete  wave  transforms 
(DWT),  transformation  plots  were  created  as  shown  in  Figures  5  and  6.  In  Fig¬ 
ure  5,  which  is  a  2-D  frequency  versus  time-color  plot,  the  black-to-red  back¬ 
ground  indicates  normal  noise  levels.  While  the  energy  in  the  raw  signal 
increases,  the  color  changes  from  the  red  to  a  hotter  color  much  like  steel  changes 
color  as  it  is  heated.  Yellow/white  indicates  higher  energy  levels.  Figure  6  is  a 
3-D  frequency  versus  time-color  plot  similar  to  the  previous  figure,  but  it  is  shown 


‘  The  acceleration  of  gravity,  32  ft/sec^ 
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Figure  5.  Wavelet  power  spectrums  versus  time 

from  another  view  angle.  The  yellow  vertical  rises  above  the  reddish  background 
indicate  rough  areas  of  the  paved  road. 

Similarly,  the  MATLAB  toolbox  was  used  to  evaluate  the  raw  data  from  the 
gravel  road  used  as  a  bypass  of  the  bridge.  There  results  are  presented  in  Fig¬ 
ures  7  and  8  for  the  data  collected  on  the  left  front  axle  of  the  LSV.  Comparison 
with  the  paved  road  test  data  clearly  indicates  the  differences  in  energy  levels 
while  riding  over  the  gravel  road  versus  the  paved  road.  Note  the  difference  in  the 
color  variations  between  the  two  plots. 


Mound  Test  Track  Data  Results:  Wavelet 
Analysis,  3/17/00 

These  are  the  results  of  data  taken  Friday,  March  17,  2000,  at  the  Mound,  LA, 
test  track.  These  data  include  the  segment  of  road  turning  off  of  Highway  80 
down  past  the  airport  to  the  bridge  detour.  The  section  of  road  was  paved  with 
asphalt  and  relatively  smooth. 
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Figure  6.  3-D  wavelet  power  spectrum  versus  time 

There  are  seven  accelerometers  located  on  the  vehicle  which  are  labeled  as 
follows:  Chi  =  Front  Left  Axle,  Ch2  =  Front  Left  Frame,  Ch3  =  Rear  Left  Axle, 
Ch4  =  Rear  Left  Frame,  Ch5  =  Front  Right  Axle,  Ch6  =  Front  Right  Frame, 

Ch7  =  Under  the  Driver’s  Seat. 

All  segments  of  the  track  were  run  at  a  near  constant  speed  of  30  mph  with 
exception  of  the  bridge  detour  area,  which  was  traversed  at  approximately 
12  mph.  All  of  the  accelerometer  outputs  are  in  engineering  units  calibrated  to  g 
forces  (g’s).  The  following  plots  are  transformations  on  the  raw  data  using 
wavelet  analysis  called  DWT. 

Beginning  with  the  first  channel,  front  left  axle,  there  are  two  plots  of  data 
associated  with  this  channel.  Figure  7  is  a  2-D  frequency  versus  time-color  plot, 
and  Figure  8  is  a  3-D  frequency  versus  time-color  plot.  It  is  very  much  the  same 
as  Figure  7  but  is  shown  at  a  different  angle.  One  can  look  at  this  figure  and 
determine  rough  spots  in  the  road. 


16 


Chapter  4  Wavelet  Analysis  for  the  Ride  Experiments 


Figure  7.  Channel  1 ,  front  left  axle,  2-D  plot 


Mound  Test  Track  Data  Results:  Wavelet 
Analysis,  3/24/00 

Bridge  Detour.  These  are  the  results  of  data  taken  Friday,  March  24,  2000, 
at  the  Mound,  LA,  test  track.  These  data  include  the  brief  segment  around  the 
bridge  detour. 

There  are  seven  accelerometers  located  on  the  vehicle  which  are  labeled  as 
follows:  Chi  =  Front  Left  Axle,  Ch2  =  Front  Left  Frame,  Ch3  =  Rear  Left  Axle, 
Ch4  =  Rear  Left  Frame,  Ch5  =  Front  Right  Axle,  Ch6  =  Front  Right  Frame, 

Ch7  =  Under  the  Drivers  Seat. 

This  segment  of  the  track  was  run  at  a  near  constant  speed  of  30  mph  with 
exception  occurring  at  the  bridge  detour  area,  which  was  traversed  at  approx¬ 
imately  12  mph.  All  of  the  accelerometer  outputs  are  in  engineering  units 
calibrated  to  g  forces  (g’s).  The  following  plots  are  transformations  on  the  raw 
data  using  wavelet  analysis  called  DWT  (discrete  wave  transform). 
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Figure  8.  Channel  1 ,  front  left  axle,  3-D  plot 


Ride  and  shock  assessment:  LeTourneau  test  track 

1  Apr  00  -  Introduction  and  test  description.  The  objective  of  this  exper  ¬ 
iment  was  to  correlate  historical  data  collected  with  the  U.S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  (currently  known  as  ERDC)  Ride  and 
Shock  Meter  and  the  high-speed  data  collection  system  that  currently  resided  on 
the  vehicle.  The  WES  Ride  and  Shock  Meter  measures  human  body  response, 
whereas  the  VIPER  system  on  the  LSV  measures  road  surface  conditions  by 
measuring  vehicular  response.  By  mounting  accelerometers  in  and  around  the 
vehicle  in  nontraditional  locations,  this  analysis  will  attempt  to  narrow  the  scope 
of  the  problem  and  lead  to  a  viable  method  for  characterizing  varying  road 
conditions. 

Seventeen  data  files  were  collected  with  the  high-speed  data  acquisition  sys¬ 
tem,  and,  as  a  comparison,  sixteen  test  sets  were  collected  to  measure  absorbed 
power  with  the  WES  Ride  and  Shock  Meter.  The  WES  Ride  and  Shock  Meter 
measures  the  absorbed  power  that  is  subjected  to  the  driver  through  means  of  an 
accelerometer  mounted  on  the  driver’s  seat  just  below  the  groin  area.  The  high¬ 
speed  system  uses  seven  accelerometers  mounted  in  various  locations  on  the 
vehicle  to  measure  acceleration.  The  accelerometer  inputs  are  recorded  through  a 
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calibrated  amplifier  system  and  the  recorded  engineering  units  are  in  g’s.  Their 
locations  are  as  follows: 

a.  Channel  1  -  Front  Left  Axle. 

b.  Channel  2  -  Front  Left  Frame. 

c.  Channel  3  -  Rear  Left  Axle. 

d.  Channel  4  -  Rear  Left  Frame. 

e.  Channel  5  -  Front  Right  Axle. 

/  Channel  6  -  Front  Right  Frame. 

g.  Channel  7  -  Under  the  Driver’s  Seat. 

h.  Channel  8  -  Speed. 

Several  courses  were  traversed  at  different  speeds.  Namely,  Ride  Course  1 
(north  and  south  at  10,  15,  and  20  mph).  Ride  Course  4  (north  and  south  at  5,  15, 
and  20  mph).  Ride  Course  6  (north  and  south  at  10  mph),  and  the  gravel  road 
were  traversed  north  and  south  at  20  mph. 

Table  1  lists  the  raw  data  files  and  their  correspondence  to  the  event  sequence 
along  with  how  the  files  are  related  as  to  distance,  direction,  and  speed.  The 
speed  at  which  these  lanes  were  actually  run  represent  a  true  and  steady  speed. 
The  slower  speeds  were  easier  to  maintain  than  the  higher  speeds. 


Table  1 

Experiment  Program 

File  Name  &  # 

Ride  Meter# 

Ride  Course 

North  /  South 

0-500/500-0 

Speed 

Letourne.prn 

Na 

4 

S 

0-500 

10 

LetourneO.pm 

Na 

4 

N 

500-0 

20 

Letournel.prn 

1 

1 

N 

500-0 

10 

Letourne2.prn 

2 

1 

S 

0-500 

10 

LetourneS.pm 

3 

1 

N 

500-0 

15 

Letourne4.pm 

4 

1 

S 

0-500 

15 

LetourneS.pm 

5 

1 

N 

500-0 

20 

Letourne6.pm 

6 

1 

S 

0-500 

20 

LetourneT.pm 

7 

4 

S 

0-500 

5 

LetourneS.pm 

8 

4 

N 

500-0 

5 

Letourne9.prn 

9 

4 

S 

0-500 

15 

LetournelO.prn 

10 

4 

N 

500-0 

15 

LetourneH.pm 

11 

4 

S 

0-500 

20 

Letourne.12.prn 

12 

4 

N 

500-0 

20 

Letourne13.prn 

13 

6 

N 

500-0 

10 

Letourne14.prn 

14 

6 

S 

0-500 

10 

Letourne.15prn 

!  Na 

6 

S 

0-500 

10 

Letourne.16prn 

15 

3  (Road) 

N 

500-0 

20 

Letourne.17prn 

16 

3  (Road) 

S 

0-500 

20 
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Results  from  ride  meter 


The  following  figures,  Figures  9  though  12,  are  typical  plots  generated  using 
the  WES  Ride  and  Shock  Meter.  The  plots  show  absorbed  power  as  a  function  of 
speed. 


Figure  9.  Absorbed  power  plot  for  Course  1 


ABSORBED  POWER  COURSE  4. 


SPEED  (MPH) 


Figure  10.  Absorbed  power  plot  for  Course  4 
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Figure  1 1 .  Absorbed  power  plot  for  Course  6 


Figure  12.  Absorbed  power  plot  for  Course  3  (Road) 


Wavelet  Analysis  of  the  Ride  Courses 

At  the  LeToumeau  ride  test  courses,  a  sufficient  number  of  tests  at  varying 
speeds  were  conducted  to  determine  the  6-watt  absorbed  power  in  the  vertical 
direction  (  Af)., )  versus  surface  roughness  (rms,  in.)  relationship  for  the  LSV, 
using  ride  courses  1,  3, 4,  and  6.  Their  data  are  presented  in  Table  2.  The  raw 
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Table  2 

Raw  Data  of  Taken  with  the  WES  Ride  and  Shock  Meter 

Test# 

Time 

Speed,  sec 

Watts 

Iso  Rms 

Peak  Acc 

Pips 

1 

31.88 

10.69 

0.53 

0.06 

0.4 

500 

2 

29.75 

11.45 

0.59 

0.06 

0.5 

500 

3 

19.74 

17.26 

1.6 

0.11 

0.68 

500 

4 

21.77 

15.65 

1.25 

0.1 

0.63 

500 

5 

15.18 

22.45 

2.58 

0.14 

0.73 

500 

6 

16.3 

20.91 

2.09 

0.12 

0.63 

500 

7 

46.34 

7.35 

0.44 

0.05 

0.5 

500 

8 

41.55 

8.2 

0.46 

0.05 

0.45 

500 

9 

20.56 

16.58 

10.4 

0.3 

1.75 

500 

10 

21.2 

16.08 

8.57 

0.28 

1.51 

500 

11 

17.08 

19.95 

14.18 

0.35 

1.89 

500 

12 

17.55 

19.42 

10.82 

0.3 

1.42 

500 

13 

29.14 

11.69 

2.83 

0.17 

0.97 

500 

14 

30.23 

11.27 

2.06 

0.15 

0.75 

500 

15 

16.73 

20.37 

0.49 

0.06 

0.44 

500 

16 

17.41 

19.58 

0.4 

0.05 

0.42 

500 

frequency  data  collected  for  the  ride  tests  were  next  subjected  to  wavelet  analysis 
as  per  the  Mound  road  data.  By  the  process  of  trail-and-error  using  the  MATLAB 
Toolbox,  those  wavelet  families  which  appear  to  be  most  applicable  are 
Daubechies  (1992)  and  Meyer  and  Ryan  (1993).  The  elevation  profile  analysis  of 
the  first  portion  of  the  500-ft  LeToumeau  Ride  Course  1  is  shown  in  Appendix  A, 
Figure  Al. 

The  elevation  profile  shown  in  Appendix  A  is  of  the  left  track  of  LeToumeau 
Ride  Course  1.  The  rms  of  this  track  of  the  course  is  0.92.  Figure  A2,  Appen¬ 
dix  A,  shows  the  detrended  version  of  the  elevation  profile. 

The  frequency  components  in  the  ride  course  themselves  as  well  as  those  of 
the  vehicle’s  response  can  be  subjected  to  wavelet  frequency  analysis.  The  fre¬ 
quency  analysis  of  the  first  part  of  the  Letoumeau  Ride  Course  1  is  shown  in 
Appendix  A.  The  graphs,  Figure  Al  through  A6,  Appendix  A,  show  a 
Daubechies  (1992)  wavelet  analysis  of  the  detrended  ride  course  elevation  den¬ 
sity.  Both  the  fixll  300  ft  and  separately  the  first  120  ft  of  the  LeToumeau  Ride 
Course  1  (from  the  data  in  Figure  A2,  Appendix  A)  are  analyzed.  A  Daubechies 
multilevel  decomposition  was  completed  using  the  second  Daubechies  wavelet 
family.  The  top  left  graph  in  Figure  A5  is  the  level  3  approximation  to  the 
spectral  density.  The  remaining  graphs  in  Figure  A5,  display  the  next  three  levels 
of  detail  in  the  wavelet  analysis.  The  ratios  of  the  magnitudes  of  the  values  on  the 
y-axis  of  the  details  versus  the  approximation  indicate  that  it  is  hard  to  charac¬ 
terize  the  ride  course  by  this  type  of  approximation.*  This  result  supports  the 


'  Methodology  for  the  wavelet  approximations  is  that  in  the  MATLAB  toolbox  book  by  Misiti 
etal.(1996). 
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conclusion  that  the  root  mean  squared  value  of  the  course  distinguishing  the 
courses  is  adequate  to  characterize  the  courses  for  experimental  analysis  of  the 
vibrational  response  of  the  vehicle.  The  wavelet  analysis  did  not  show  that  a 
significant  part  of  the  frequency  power  spectrum  of  the  course  was  left  in  the 
details  of  the  first  three  levels  of  approximation  to  the  signal.  Thus,  one  would 
expect  a  Fourier  transform  or  wavelet  transform  to  not  be  affected  by  the  response 
of  the  vehicle’s  vibrations  to  this  part  of  the  input.  And,  this  supports  the  notion 
that  the  root  mean  square  statistic  is  adequate  to  characterize  the  course  for 
experimental  analysis. 

Figures  13  and  14  show  the  results  of  harmonic  wavelet  analysis  of  the  LSV’s 
test  response  to  Letoumeau  Ride  Course  1. 


Figure  13.  Channel  7,  2-D  plot  for  acceleration  as  frequency  versus  time 

This  harmonic  wavelet  analysis  plot  indicates  that  most  of  the  harmonic 
wavelet  energy  of  the  response  of  the  vehicle  to  the  ride  course  is  concentrated  in 
the  low-frequency  area.  With  respect  to  the  spring  constants  in  the  vehicle  design, 
the  natural  frequency  of  the  test  vehicle  frame  is  expected  to  be  less  that  5HZ. 

This  natural  frequency  can  be  determined  experimentally  using  a  drop  test  and 
measuring  the  time-decay  constant  of  the  plot  of  how  the  springing  motion  damps. 
These  plots  suggest  that  it  is  possible  to  characterize  the  vehicle’s  response  to  the 
course  using  a  few  low-order  wavelet  coefficients.  To  do  this,  it  is  necessary  to 
use  a  frequency  ordered  wavelet  decomposition  tree  of  the  signal. 
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Figure  14.  Channel  7,  3-D  plot  for  acceleration  as  frequency  versus  time 

Figure  15  shows  the  Meyer  wavelet  decomposition  tree  and  frequency  ordered 
wavelet  coefficient  values  for  the  first  500  data  points  of  the  Letoumeau  Ride 
Course  1  data.  The  x-axis  of  the  box  (Figure  15,  in  box  at  lower  right)  is  time  and 
the  y-axis  is  frequency.  The  wavelet  decomposition  tree  has  been  cut  off  at  the 
second  level  of  decomposition.  The  number  of  boxes  on  the  vertical  y-axis 
corresponds  to  the  number  of  frequency  decomposition  nodes  on  the  lowest  level 
in  the  wavelet  decomposition  tree.  The  results  show  that  most  of  the  wavelet 
energy  is  concentrated  in  the  first  node  of  the  wavelet  packet  decomposition  tree. 

The  analysis  indicates  that  this  test  vehicle’s  vibrational  response  to  the  ride 
course  can  be  characterized  by  the  value  of  one  Meyer  wavelet  coefficient. 

Figure  16  shows  the  automatic  de -noising  of  the  Letoumeau  ride  course  ride 
meter  data  using  a  data  compression  method  called  soft  wavelet  thresholding.  The 
graph  in  upper  right  box  (Figure  16)  in  the  yellow  color  is  the  reconstmcted 
de-noised  signal.  The  de-noising  threshold  applied  was  applied  to  cutoff  wavelet 
coefficients  whose  sorted  absolute  value  was  below  0.77  of  the  maximum 
coefficient.  This  threshold  was  computed  from  a  formula  given  in  the  book  by 
Misiti  et  al.  (1999)  which  involved  a  logarithm  of  the  number  of  data  points.  As 
mentioned  above,  in  this  case,  most  of  the  wavelet  packet  energy  is  concentrated 
in  the  first  node  of  the  wavelet  packet  decomposition  tree.  Thus,  thresholding 
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Cut  Tree  Analyzed  Signal ;  length  =  3200 


100  200  300  400  500 


Colored  Coefficients  for  Terminal  Nodes 


frequency  ordered  coefficients 


.  .  '  ■  - 

Scale  of  Colors  from  Min  to  Max 


Figure  15.  Meyer  Wavelet  Decomposition  Tree 

seems  appropriate  here.  The  x-axis  of  the  box  (Figure  16,  box  at  lower  right)  is 
time  and  the  y-axis  is  frequency.  Higher-order  frequency  components  have  been 
dropped  from  the  signal. 

For  each  of  the  3-D  plots,  such  as  Figure  6  or  8,  it  is  possible  to  calculate  the 
value  of  the  area  under  the  plotted  surface.  This  value  is  in  G’s.  The  amplitude, 
or  values  plotted  against  the  z-axis,  is  also  in  G’s.  The  surface  is  actually  a 
function  called  A(t,f)  where  A  is  the  amplitude  of  acceleration,  t  is  time,  and / 
is  frequency.  Integrating  acceleration  with  respect  to  time  yields  velocity.  Inte¬ 
grating  A(t,f)  ,  with  respect  to  time  and  frequency,  yields  the  original  unit, 
because  the  units  of  frequency  are  Hz  =  1/ 1 . 


Chapter  4  Wavelet  Analysis  for  the  Ride  Experiments 


25 


Figure  16.  Ride  course  threshold  coefficients 


V^ffA(t,f)dtc//  =  Vtl/t  (10) 

Using  these  numbers  and  two  selected  channels,  specifically  channels  3  and  5, 
some  very  interesting  features  appear.  The  total  volume  numbers  from  Courses  1 , 
4,  and  3  (the  gravel  road)  were  calculated.  The  total  volume  decreased  as  speed 
increased.  This  is  to  be  expected  because  the  length  of  time  through  the  test  track 
decreases  as  speed  increases.  Dividing  the  total  volume  g’s  by  the  time  in  the 
track,  or  the  length  of  the  record,  the  g’s  were  for  the  most  part  constant,  regard¬ 
less  of  speed.  The  actual  tabulated  data  are  given  in  Table  3,  and  the  plots  are 
shown  in  Figures  17  and  18. 
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Table  3 

G/sec  Calculation  Table 

Course  1 

Ch3 

Ch5 

Rms=0.77 

Rear  Left  Axle 

Front  Right  Axle 

Time 

Total  g’s 

Total  g’s 

Total  g's 

Total  g's 

sec 

Speed 

N 

S 

N 

S 

N 

S 

10  mph 

212.13 

205.96 

162.52 

165.59 

32.00 

29.50 

15  mph 

155.42 

172.04 

120.12 

128.48 

19.50 

22.00 

20  mph 

111.81 

130.34 

91.22 

101.33 

15.00 

16.50 

g’s/sec 

g's/sec 

g's/sec 

g's/sec 

10  mph 

6.63 

6.98 

5.08 

5.61 

15  mph 

7.97 

7.82 

6.16 

5.84 

20  mph 

7.45 

7.90 

6.08 

6.14 

Avg 

7.35 

7.57 

5.77 

5.86 

Course  4 

Ch3 

Ch5 

Rms=1.16 

Rear  Left  Axle 

Front  Right  Axle 

Time 

Total  g's 

Total  g's 

Total  g’s 

Total  g's 

sec 

Speed 

N 

S 

N 

S 

N 

S 

10  mph 

320.18 

261.26 

269.36 

209.45 

46.50 

41.50 

15  mph 

208.30 

238.38 

156.35 

177.38 

20.50 

21.00 

20  mph 

528.35 

550.25 

454.02 

494.60 

60.50 

45.00 

g's/sec 

g's/sec 

g’s/sec 

g's/sec 

10  mph 

6.89 

6.30 

5.79 

5.05 

15  mph 

10.16 

11.35 

7.63 

8.45 

20  mph 

8.73 

12.23 

7.50 

10.99 

Avg 

8.59 

9.96 

6.97 

8.16 

Course  3 

Ch3 

Ch5 

Rms=0.5 

Rear  Left  Axle 

Front  Right  Axle 

Total  g's 

Total  g's 

Total  g's 

Total  g's 

Time 

Speed 

N 

S 

N 

S 

N 

S 

20  mph 

86.13 

82.62 

64.64 

66.35 

16.5 

17.5 

g's/sec 

5.22 

4.72 

3.92 

3.79 

Rear  Left  Axle 

Front  Right  Axle 

RMS 

N 

S 

N 

S 

avg 

avg 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.50 

5.22 

4.72 

3.92 

3.79 

4.97 

3.85 

0.77 

7.35 

7.57 

5.77 

5.86 

7.46 

5.82 

1.16 

8.59 

9.96 

6.97 

8.16 

9.28 

7.57 
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The  values  of  interest  here  are  g’s/sec.  These  values  are  plotted  in  the  fol¬ 
lowing  figures:  Figure  17,  Course  1,  in  which  the  g’s/sec  are  plotted  as  a  function 
of  speed,  as  well  as  for  Figure  18,  Course  4.  The  results  for  Course  1  show  that 
the  g’s/sec  are  fairly  constant  with  the  differences  being  channel  to  channel.  RLA 
is  Right  Left  Axle  and  FRA  is  Front  Right  Axle.  N  and  S  are  North  and  South. 
Course  1  appeared  to  be  uniform  regardless  of  which  direction  it  was  traversed. 
This  did  not  appear  to  be  the  case  with  Course  4. 

The  reason  Course  3  was  included  in  this  analysis,  even  though  it  presents 
only  one  point,  was  in  order  to  use  a  lower  rms  value.  Averaging  the  g’s/sec 
values  so  that  one  value  is  representative  for  the  speeds  given  in  our  test  case 
yields  the  Avg  values  shown  in  Table  3.  Taking  these  average  values  and  pre¬ 
senting  these  values  given  in  the  lower  portion  of  Table  3,  plotting  rms  with 
respect  to  g’s/sec  results  in  Figure  19. 


This  particular  plot  shows  a  second-order  polynomial  curve  fit  through  the 
data.  But  keep  in  mind  that  it  is  only  valid  for  limits  of  the  data  set  presented.  A 
broader  scope  could  be  obtained  with  different  courses  and  at  higher  speeds  with 
more  experienced  drivers. 
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(11) 


rms  =  0.0077 (g's/sQcf  +  0.0928(g'5/sec)  +  0.004 

For  a  single  channel  only. 

Given  a  broad  database  of  varying  vehicles’  ride  performance  data,  one  could 
make  a  sound  judgment  on  the  road  condition’s  impact  to  that  vehicle.  The  road 
could  then  be  categorized  and  speed  restrictions  placed  for  that  section  of  road. 
Figure  20  shows  the  absorbed  power  data  collected  from  the  ride  meters  as  a 
function  of  rms,  which  is  a  family  of  curves  for  this  vehicle.  A  5-ton  truck  will 
obviously  have  a  different  family  of  curves,  but  the  constant  that  remains  is  the 
indirectly  measured  rms  value. 


Ride  Curves 


— I— Seriesi 
— Series2 
SeriesS 
— Series4 
— SeriesS 
— Series6 


Figure  20.  Absorbed  power  curves  for  Camouflage  Tactical  Utility  Vehicle 
(CTUV)  Suburban  vehicle 
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Conclusions 


A  Meyer  wavelet  decomposition  tree  and  frequency-ordered  wavelet  coeffi¬ 
cient  values  were  computed  for  the  first  500  data  points  of  the  LeToumeau  Ride 
Course  1  data.  The  wavelet  decomposition  tree  was  cut  off  at  the  second  level  of 
decomposition.  The  results  showed  that  most  of  the  wavelet  energy  is 
concentrated  in  the  first  node  of  the  wavelet  packet  decomposition  tree.  Thus,  the 
analysis  indicates  that  for  this  ride  course,  this  test  vehicle’s  vibrational  response 
to  the  ride  course  can  be  characterized  by  the  value  of  one  Meyer  wavelet 
coefficient. 

The  first  part  of  the  LeToumeau  Ride  Course  1  was  compared  in  terms  of  its 
Fourier  power  spectral  transform  and  wavelet  power  spectral  transform.  The 
analysis  indicated  no  need  to  use  wavelet  analysis  to  characterize  the  ride  course 
itself  This  result  supports  the  conclusion  that  the  root  mean  squared  value 
distinguishing  the  course  is  adequate  to  characterize  the  course  for  experimental 
analysis  of  the  vibrational  response  of  the  vehicle.  The  wavelet  analysis  did  not 
show  that  significant  parts  of  the  frequency  power  spectmm  of  the  course  were 
left  in  the  details  of  the  first  3  levels  of  approximation  to  the  signal.  Thus,  one 
would  expect  a  Fourier  transform  or  wavelet  transform  not  to  be  affected  by  the 
response  of  the  vehicle's  vibrations  to  this  part  of  the  input.  And,  this  supports  the 
notion  that  for  the  off-road  or  cross-country  ride  courses  the  root  mean  square  sta¬ 
tistic  is  adequate  to  characterize  the  course  for  these  experiments. 

The  harmonic  wavelet  analysis  plot  of  the  test  vehicle’s  vibrational  response 
to  the  LeToumeau  Ride  Course  1  indicates  that  most  of  the  harmonic  wavelet 
energy  of  the  response  of  the  vehicle  to  the  ride  course  is  concentrated  in  the  low- 
frequency  area.  With  respect  to  the  spring  constants  in  the  vehicle  design,  the 

natural  frequency  of  the  test  vehicle’s  frame  is  expected  to  be  less  that  5HZ.  This  / 

natural  frequency  can  be  determined  experimentally  using  a  drop  test  and 
measuring  the  time  decay  constant  of  the  plot  of  how  the  springing  motion  damps. 

These  plots  suggest  that  it  is  possible  to  characterize  the  vehicle’s  response  to  the 
course  using  few  low-order  wavelet  coefficients.  To  do  this,  it  is  necessary  to  use 
a  frequency-ordered  wavelet  decomposition  tree  of  the  signal. 

The  Mound  Test  Road  results  indicate  that  wavelets  can  be  used  to  analyze 
the  vehicle  accelerometer  data  plotted  with  respect  to  time.  A  3-D  surface  plot  of 
the  data  can  also  be  used  as  a  measure  of  the  total  absorbed  power  by  the  vehicle 
over  the  ride  course.  This  plot  provides  an  alternate  means  of  computing  this  other 
than  using  the  ride  meter,  which  has  been  in  use  for  some  years  now.  Finally,  it 
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should  be  mentioned  that  road  test  tracks  will  not  necessarily  have  simply 
designed  power  spectrums  like  the  ERDC  off-road  courses.  Thus,  further 
research  needs  to  be  conducted  to  evaluate  the  usefulness  of  using  wavelets  to 
characterize  road  conditions. 
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Appendix  A 

Wavelet  Analyses  of  the 
LeTourneau  Ride  Courses 


The  graphs,  Figures  A1  through  A6,  in  this  appendix  show  a  Daubechies 
(1992)^  wavelet  analysis  of  the  detrended  ride  course  elevation  density.  The  first 
120  ft  and  300  ft  of  the  LeTourneau  1  Ride  Course  (data  in  Figure  Al)  were  ana¬ 
lyzed  separately.  Figure  A2  shows  the  detrended  values  of  the  elevations  in  the 
same  course.  Figure  A3  shows  the  power  spectral  density  computed  from  the 
detrended  elevations.  Figure  A4  shows  a  Daubechies  multilevel  decomposition 
using  the  second  Daubechies  wavelet  family.  Top  left  graph  of  Figure  A4  is  the 
level  3  approximation  to  the  spectral  density.  The  remaining  graphs  display  the 
next  three  levels  of  detail  in  the  wavelet  analysis.  The  ratios  of  the  magnitudes  of 
the  values  on  the  y-axis  of  the  details  versus  the  approximation  indicate  that  it  is 
hard  to  characterize  the  ride  course  by  this  type  of  approximation.^ 

The  term  rms  (root  mean  squared  elevation,  inches)  is  used  by  the  Mobility 
Systems  Branch,  Geotechnical  and  Structures  Laboratory  (GSL),  U.S.  Army 
Engineer  Research  and  Development  Center  (ERDC),  Vicksburg,  MS,  site,  in  the 
characterization  of  the  surface  roughness  of  terrain.  It  is  determined  by  first 
detrending  (subtracting  the  linearly  increasing  or  decreasing  component)  the 
surface  elevation  measurements  taken  at  1-ft  intervals  in  the  terrain  profile  and 
then  computing  the  ordinary  square  root  of  the  variances  of  the  measurements 
from  the  detrended  value.  The  paper  by  Harrell  (1966)  gives  a  fuller  discussion  of 
how  rms  is  used  in  vehicle  tests  to  determine  ride-limiting  speed.  Ride-limiting 
speed  represents  the  speed  at  which  vehicle  vibrations  at  the  driver’s  seat  reach 
absorbed  power  limits  of  6  watts  in  the  vertical  direction. 

Figure  A3  depicts  the  power  spectral  density  of  the  profile  elevation  in  the 
first  120  ft  of  the  LeTourneau  Ride  Course  1  (left  track).  The  graphs  plotted  hori¬ 
zontally  up  to  the  limits  of  the  nyquist  sampling  frequency  (50  cycles/second)  of 
the  instrumentation  on  the  test  vehicle.  The  values  were  calculated  by  taking 
discrete  Fourier  transforms  of  the  data. 


^  References  are  listed  at  the  end  of  the  main  text  of  this  report. 

2 

Methodology  for  the  wavelet  approximations  is  that  in  the  MATLAB  toolbox  book  by  Misiti  et  al 
(1996). 
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A1 


Figure  A1 .  Elevation  profile  of  LeTourneau  Ride  Course 


Figure  A2.  Detrended  elevation  profile 
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Figure  A3.  Power  spectral  density  of  the  ride  course  elevation  profile 
LeTourneau  1  (left  track) 


Figure  A4.  Wavelet  analysis  of  ride  course  1  (0  to  120  ft) 
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Wavelet  Approximation  of  Detrended  Ride  Course 
Elevations 

The  graphs  shown  in  Figure  A4  and  A5  show  a  Daubechies  (1992)  wavelet 
analysis  of  the  detrended  ride  course  elevation  density.  Again,  the  first  120  ft  of 
the  LeToumeau  Ride  Course  1  (Figure  A4)  and  the  entire  300  feet  (Figure  A5) 
were  analyzed.  A  Daubechies  multilevel  decomposition  was  accomplished  using 
the  second  Daubechies  wavelet  family.  The  top  left  graph  is  the  level  3 
approximation  to  the  spectral  density.  The  remaining  graphs  display  the  next 
three  levels  of  detail  in  the  wavelet  analysis  (Dl,  D2,  and  D3).  The  ratios  of  the 
magnitudes  of  the  values  on  the  y-axis  of  the  details  versus  the  approximation 
indicate  that  it  is  difficult  to  characterize  the  ride  course  using  detrended 
elevations  with  this  type  of  approximation.' 


Approximation  A3  Detail  D1 


Figure  A5.  Wavelet  analysis  of  ride  course  1  (0  to  300  ft) 


’  Methodology  for  the  wavelet  approximations  is  that  in  the  MATLAB  toolbox  book  by 
Misti  et  al.  (1996),  listed  in  references  following  main  text. 


A4 
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Wavelet  Approximation  of  Ride  Course  Power 
Spectral  Density  Curve 


Figure  A6  shows  the  same  analysis  as  the  previous  figures  except,  in  this  case, 
the  same  wavelet  families  were  used  to  approximate  the  power  spectral  density 
graphs. 

Here  again,  comparing  the  magnitude  of  the  detail  coefficients  to  the  level  3 
approximation  coefficients,  the  approximation  does  not  adequately  capture  the 
data  to  the  detail  desired. 
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Decomposition  at  ievel  3  :  s  =  a3  +  d3  +  d2  +  d1  . 
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